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Abstract 

The three-state agent-based 2D model of financial markets in the version 
proposed by Giulia Iori in 2002 has been herein extended. We have intro- 
duced the increase of herding behaviour by modelling the altering trust of 
an agent in his nearest neighbours. The trust increases if the neighbour has 
foreseen the price change correctly and the trust decreases in the opposite 
case. Our version only slightly increases the number of parameters present 
in the Iori model. This version well reproduces the main stylized facts ob- 
served on financial markets. That is, it reproduces log-returns clustering, 
fat-tail log-returns distribution and power-law decay in time of the volatility 
autocorrelation function. 



PACS numbers: 89.65. Gh, 02.50.Ey, 89.75.Fb, 45.70.Vn. 



1 Introduction 

The modern study of financial markets has discovered a huge number of phenomena 
violating the Brownian stochastic dynamics of financial markets. Three major styl- 
ized facts observed in the market are exceptionally intriguing: volatility clustering, 
fat-tail log-return distribution and power-law decay in time of the volatility auto- 
correlation function. The promising approach to reproduce these phenomena runs 
*j. lipski@student.uw.edu. pi 



through heteroagent-based models where decisions of individual agents somehow 
determine the price dynamics. Various versions of such models were thoroughly 
discussed in literature [1-7]. Our version well describes the above mentioned styl- 
ized facts by introducing a characteristic emotion to the Iori model. That is, the 
trust in the foreseeing market agents. 

2 Outline of the Iori model 

The Giulia Iori model [1] consists of n spins or agents placed in 2D square lattice 
sites, assuming one of the following values of spin states: either +1 as buy, —1 as 
sell or as stay passive. Each agent, % = 1, 2, n, is under the influence of a local 
field Yi which is the sum of the forces exerted by the nearest neighbours and by 
random noise. 

Thresholds play a significant role in the Iori model. Thanks to these thresholds, 
the temporal spin <7j is defined as follows: 



where is a time-dependent positive threshold considered below and Yi(t) is a 
time-dependent local field: 



here 3^ is a time-dependent force exerted by agent j on his neighbouring agent 
i and T]i(t) is a time-dependent white noise describing the agent's emotion. This 
force is defined in Section 13.11 

The threshold reflects the symmetry between the probability of buying 
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and selling. The initial threshold £i(0) is randomly drawn from the Gaussian 
distribution N(0, 1). Next, it is adjusted in successive time steps after each decision 
round, according to the rule: 

&(* + *) - P ® t>0 (3) 

where P(t) is a market price of a share at time t. For the negative time (i.e. for 
history before the begining of the simulation) the price values P(t < 0) are drawn 
from a uniform distribution. 

In Equation ^ it was tacitly assumed: (i) the proportionality of the threshold 
to the transaction cost and (ii) proportionality of this transaction cost to the share 
prices. Iori proved that if the threshold is either zero or constant in time, the 
stylized facts cannot be properly reproduced. 

It is assumed that agents mutually exchange information in a consultation 
round and they have the opportunity to change their opinion (or spin state) once 
on average. Hence, the local field can alter. The consultation round is the analogy 
to the onr Monte Carlo step per spin in Monte Carlo simulations of magnetism. 
The agents begin trading only after the field is relaxed. It is empirically proven 
within the Iori model that then the most of spins become constant, which is an 
analogy to thermalisation of magnetisation in physics. 

Subsequently, one calculates the total supply S(t) (the total number of negative 
spins) and demand D[t) (the total number of positive spins). Thus, the market 
maker can determine the current market price according to the following rule: 

p(t+i) (D(t)y® 

P{t) \S{t)) ' u 

where market activity 

D(t) + S(t) 

Kit) = a (5) 

n 



is a slowly-varying function of time and coefficient a was assumed as a sufficiently 
small calibration parameter. Equation Q describes an asymmetric reaction of 
the market maker to imbalance between supply and demand. The intensity of 
this reaction (measured by the value of log-return) depends linearly on the market 
activity K,(t). Equation Q is consistent with observed positive correlation between 
absolute log-returns and trading volume (defined, as usual, by quantity related to 
min(D(t), S(t))). The relation between demand, supply and the price change 
is differently considered in various agent-based models [6 J (and refs. therein). 
However, the models have to be consistent with the classic economic law of demand 
and supply. 

As usual, the log-return r over a time period t — t' is defined as a natural 
logarithm of ratio of the corresponding prices P(t) and P(t — £'): 

In fact, mainly this definition is exploited in our paper. 

3 Our extension: trust in the foreseeing neighbours 
3.1 Definition of the impact 

The aim of our extension is to define a more realistic impact of a given agent 
on his neighbours than in the Iori model. This impact is large if the agent has 
recommended to buy the asset before a price increase or to sell before a price 
decrease. That is, the impact is large if the product of the spin value at a certain 
time in the past and the log-return from that time until now is positive. Otherwise, 
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the impact is small. Hence, the force exerted by an agent j on his neighbour % is: 

Mt) = w i3 + £ ln (^v) > r ^ 2 - ( 7 ) 

Coefficient Jjj(i = 0) is the value of the initial force allotted at the beginning of the 
simulation to the pair of investors It is either 1 with a fixed probability p or 
with complementary probability 1—p (as in Iori model) and Wij is the background 
static influence of agent j on agent i, randomly drawn from a uniform distribution. 
The sum over r' represents the dynamic part of the impact containing a kind of 
r-step memory. 

To avoid a persistent positive feedback effect (which results in directed constant 
price changes), a fundamental behaviour of agents is introduced. This fundamental 
behaviour is constrained by the fundamental price and two positive factors a(> 1) 
and b(< 1). If the market price is greater than the fundamental price multiply by 
a, an agent sells shares. Reversely, if the market price is lower by factor b, an agent 
buys shares. 

3.2 Algorithm 

Our algorithm makes possible to simulate the behaviour of agents. Single agent 
can trade only one share at a time that is, he can buy shares only if he possesses 
enough cash or sell shares only if he has at least one share. Leverage and short 
sale is beyond the scope of our model. 

The algorithm consists of the initial and proper parts. Within the initial 
part the input values of variables and required parameters are prepared. This 
part is valid for the case t — and — r < r' < —1, where r = 1, 2, . . . . For each 
value of t' (at fixed value of r) it consists of the following steps. 

1. The spin state, <7i(r') (i = 1, 2, n), of each agent is drawn from a (discrete) 
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3-point uniform distribution. 

2. The price P(r') is drawn from unit interval, (0, 1), by using a uniform con- 
tinuous distribution. 

3. The fixed values W%j (valid for any time) are drawn from interval (—2, 2) by 
using a uniform distribution. 

4. Thanks to the above given steps we can calculate Jij(t = 0) from Equation 
Q and subsequently Yi(t = 0) from Equation [2j 

5. The intial values of the thresholds ^(t = 0) are drawn from the normal 
distribution JV(0, 1). 

6. The initial values of the 5-correlated noise, r]i(t = 0), are drawn from the 
normal distribution N(0, 1). 

7. Cash and shares are granted as well to every agent as to market maker (for 



details see Section 3.3). 



Notably, each time step is divided into two rounds: the consultation and de- 
cision ones. In the consultation round only the relaxation of spins occur without 
changing any amount of cash and shares (of each agent and market maker). The 
change of cash and shares takes place only in the decision round. This is described 
in details below. 

Within the proper part, valid for t > 1, the dynamics of the system is 
simulated. In this part the dynamics of every agent is determined by the following 
algorithm. 

1. For t = 1 the spin state of each agent is calculated according to Equation 
Q, as we have already all quantities required (i.e. those for t — 0). Thus, 
for t = 1 all decisions of agents are known. 
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2. Next, agent i is drawn with probability 1/n. 

3. By using the spin values of i's nearest neighbours, we construct forces Jy(t = 
1) exerted on agent % (from his nearest neighbours js) according to Equation 
([7]). This is possible as all required quantities have been calculated one step 
earlier. 

4. Hence, the local field Yi(t = 1) is calculated according to Equation §2§. 

5. The agent's threshold = 1) is calculated from Equation (|3| as all required 
quantities have been calculated one step earlier. 

6. Finally, the spin state crj(t = 2) is calculated according to Equation Q. 

7. The steps 2-6 are repeated until the spin relaxes. This means that decisions 
of agents stabilise and the consultation round is finished. This relaxation is, 
in fact, observed in our simulation quite well. 

8. The final spin values are the final decisions of all agents. That is, these spin 
values are the output of the consultation round preceding the decision round 
in which agents are trading. 

9. According to their final decisions (taken at the end of the consultation round), 
agents put sell or buy orders (being only allowed to trade one share at a given 
time step) however, the final share prices still have to be defined. 

10. The market maker determines the price according to Equation |4} 

11. The agents trade, as long as they possess enough shares or cash. 

12. Agents who have no sufficient amount of cash or shares, are not able to trade 
with other agents. Instead, they can trade with the market maker (at prices 
established in above item) if he has sufficient amount of cash and shares. 
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However, if market maker has no sufficient cash and shares the algorithm 
stops the current decision round and goes to the subsequent time step. 

13. In a given decision round t (the number of decision round always equals the 
number of time steps) the agents demonstrate a fundamental behaviour with 
probability ir under the following conditions: 

• remainder R from dividing the time step (or number) t (t — 1, L) by 
a natural number m(= K + k) is lower than a fixed value p (or R = 
t — Ent(t/m)m < p, where Ent(x) means the entire part of x). Number 
K is a natural number fixed at the beginning of the simulation, while k 
is also a natural number but drawn from a uniform distribution. As k 
fluctuates, it protects the agent trading against periodicity (which could 
be present in the case of fixed m). Apparently, for t < K remainder R 
is always smaller than p (as it vanishes); 

• the market price of share is either higher than the fundamental price 
multiplied by factor a(> 1) (then the agent can sell a share for the 
market price) or lower than the fundamental price by factor 6(< 1) 
(then the agent can buy a share for the market price). Indeed, this is 
a fundamental behaviour of agent because he is driven by the relation 
between the market price and the fundamental one and not by the 
collective impact of his neighbours. 

3.3 Comments concerning simulations 

The values of parameters of the model have a substantial influence on the accu- 
racy of the reconstruction of stylized facts. For instance, decreasing coefficient 
a (present in Equation (|5])) by one order of magnitude (here from a = 0.01 to 
a = 0.001), increasing the number of agents n (here from n = 10 to n = 15), 
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and decreasing probability it of trading according to the fundamental price (here 
from 7r = 90% to 7r = 80%), one gets less accurate results: less visible log-returns 
clustering and less accurate log-returns distribution. 

We performed our simulations, named simulation A and simulation B, which 
vary only by values of some parameters. The parameter values of simulation A are 
as follows: 

• the number of agents n = 1024; 

• at the beginning of the simulation each agent a'priori received 100 cash units 
and the same number of shares. Similarly, the market maker received 10240 
cash and shares; 

• L = 80000 decision rounds were set and the system's memory was extended 
until t = 20 decision rounds; 

• the agents demonstrated a fundamental behaviour (see Algorithm in Section 



3.2 for details) if the remainder of dividing the number of decision round t by 
m = 275 + A; (where K = 275 and discrete variable k was drawn from domain 
k — 1, . . . , 15 by a uniform distribution) was lower than p = 20; thus they 
sell shares if their price is a = 1.5 times greater than the fundamental price 
and buy shares if their price is b = 0.667 times lower than their fundamental 
price - both with probability 7r = 70%; 

• in every decision round the fundamental value of shares rise by factor j§^j) 

• the coefficient a from Equation Q equals 0.01. 

The values of parameters driving simulation B were assumed the same as in 
simulation A, except the following ones: 

• the system's memory r = 40 decision rounds; 
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the agents demostrated fundamental behaviour (see Algorithm in Section 



3.2 for details) if the remainder of dividing the number of decision round t 
by m = 200 + k (where K = 200 and discrete variable k was drawn from 
domain 1, . . . , 10 by a uniform distribution) was lower than p = 30; then they 
traded according to the fundamental price analogously as for simulation A 
(considered above) but with larger probability 7r = 90%. 



For the purpose of comparison, we have also conducted simulations for the case 
where the agents' memory was extended only over one time i.e. over r = 1. For 
comparison, these simulations with no agents' esteem, are shown together with 
simulations A and B in Figure [TJ There the long-term rising and falling of trust in 
the foreseeing agents is somehow implicitly coded. 



4 Main results 

4.1 Log-returns clustering 

The highly non-Gaussian property of log-returns as defined by Equation (|6| of log- 
returns time clustering is being studied in this section. The results given below, 
showing log-returns clustring, have been obtained by assuming 6 decision rounds 
as a single trading day. This behaviour directly implies the clustering of volatility 
(defined, e.g., as the absolute value of log-returns or their square). Our results 
obtained in simulations A and B, are shown in Figure [I] (both middle plots), were 
compared with the closing time daily empirical data for the indicies: the Warsaw 
Stock Index (WIG) from 1991-04-16 to 2012-01-05 (top left plot) and the Standard 
& Poor's 500 index from 1968-01-01 to 2012-01-05 (top right plot) [10]. 

Apparently, the log-returns clusterings obtained in simulations for variograms 
are sufficiently distinct, although not so pronounced as for the corresponding em- 
pirical variograms. All variograms are shown in Figure [TJ Notably, no periodicity 
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is observed for the variograms despite the fact that the fundamental behaviour of 
investors has some periodic component. Our results occur not to be substantially 
affected if the system's memory is cut to the minimal range, that is to r = 1 (see 
the plot placed at the bottom of Figure [TJ. 

4.2 Power-law decay of autocorrelation function 

The (normalized) autocorrelation function over time t as a function of the time lag 
r has been calculated as follows: 

= (r(t)r(t-r))-(r) 2 ^ 
Var(r) 

where t is the decision round number, r is the time lag, (...) is a time average (cf. 
Equation (pi)), and Var(r) is the variance of the distribution of that log-return 
over time period t. 

The predictions of Equation ^ are compared with empirical data in Figure [2j 
Despite of poor agreement, some short-term relaxations are well seen, in particular, 
for the S&P 500 index (bottom plots). 

The autocorrelation function of absolute daily log-returns reveals long-term 
power-law relaxation versus time: 

R(r) oc r" 7 , 7 > 0. (9) 

The values of exponent 7 obtained from simulations well agree with the correspond- 
ing ones obtained from the empirical data both for WIG and S&P 500. We found 
the following values of exponent 7 for the empirical data namely, 7 (WIG) = 0.546 
and 7(S&P500) = 0.541 as well as for simulations 7(A) = 0.576 and 7(5) = 0.503 
(cf. Figure [3]). Exponent 7 for simulation A agrees well with the corresponding one 
for WIG data and for simulation B agrees with that for S$P 500 data. Moreover, 
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the exponent 7 for the model without long memory of agents' esteem (i.e. for 
t = 1) equals j(t = 1) = 0.504, which is not so far placed from results obtained 
in simulations A and B. The results shown in Figure | make our model promising 
for studying subsequent stylized facts. 

4.3 Fat-tail log-returns 

Histograms for real market daily log-returns exhibit fat-tails (see Figure [4]). How- 
ever, the tails are thinner than the corresponding tails of the Levy distribution. 
Apparently, better agreement between simulations and empirical data has been 
observed for the developed market described by S&P 500 than for the emerging 
Warsaw Stock Exchange. This stylized fact (i.e. the presence of fat-tails in daily 
log-returns distributions) is not reproduced if the system's memory is assumed as 
t = 1. 

5 Conclusions 

Our modifications of the Iori model using success based strategy, which introduced 
intrinsically-driven herding and emotional behaviour of investors, gave satisfactory 
results. The following stylized facts have been reproduced: 

• log-returns and hence volatility clustering, 

• fat-tail log-return distribution, 

• power law autocorrelation function decay vs. time of log-returns and also 
quick (short-term) decay of autocorrelation function of log-returns. 

It is shown that the role of a long memory of the system, that is, a growing trust 
in the foreseeing neighbours, is most pronounce in the case of reproducing fat-tail 
dostribution of log-returns. 
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We can conlude our (as we hope) quite representative results by the general 
observation that our model significantly better fits empirical data coming from 
stock markets of large size than from emerging markets. 

Since the local field YJ (present in Equation ([2])) can be interpreted as the ith 
agent's utility function, a bridge between our model and utility function theories 
has been constructed. 
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Daily log-returns clustering for real market [TO] and for all our sim- 
ulations. Detailed descriptions of the plots are in their titles and 

legends [15] 

Autocorrelation function for daily log-returns for real markets [10] 
and for our simulations A and B vs. time. The autocorrelation 
function is less oscillating in our simulations than for the empirical 
data. However, in both cases its short-term relaxation is well seen. 
Detailed descriptions of the plots are in their titles and legends. . . [16] 
The empirical (pluses) and simulational (crosses) autocorrelation 
functions of absolute log-returns vs. time. Autocorrelation func- 
tions obtained from simulations are consistent with the correspond- 
ing ones for market data |1U| for both indicies. For comparison, 
simulations with no agents' esteem are also shown. The power-law 
relaxations of all autocorrelation functions are well seen over several 
decades. Detailed descriptions of the plots are in their titles and 

legends [T71 

Histograms of daily log-returns both for market data [TU] and ob- 
tained from simulations in the semi-logarithmic plots. The obtained 
results show the non-Gaussian properties of the price process, which 
is satisfactorily described by our model, in particular, for the stock 
market of a large capitalisation. Detailed descriptions of the plots 
are in their titles and legends [TH] 
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Daily log-returns, WIG data since 1991 vs. simuiation A 
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Figure 1: Daily log- returns clustering for real market [10] and for all our simula- 
tions. Detailed descriptions of the plots are in their titles and legends. 
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Autocorrelation of log-returns, simulation A v. WIG data since 1991 




-0.01 1 ■ — ■ — 1 

1 10 100 1000 10000 

time lags [trading days] (log scale) 
Autocorrelation of log-returns, simulation B v. WIG data since 1991 




-0.01 1 ■ — ■ — 1 

1 10 100 1000 10000 

time lags [trading days] (log scale) 



Autocorrelation of log-returns, simulation A v. S&P50Q data since 196£ 




-0.003 1 ■ — ■ — 1 

1 10 100 1000 10000 

time lags [trading days] (log scale) 
Autocorrelation of log-returns, simulation B v. S&P5QQ data since 1965 




-0.003 1 ■ — ■ — 1 

1 10 100 1000 10000 

time lags [trading days] (log scale) 



Figure 2: Autocorrelation function for daily log-returns for real markets [10] and for 
our simulations A and B vs. time. The autocorrelation function is less oscillating 
in our simulations than for the empirical data. However, in both cases its short- 
term relaxation is well seen. Detailed descriptions of the plots are in their titles 
and legends. 
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Autocorrelation of volatility, simulation A v. WIG data since 1991 
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Autocorrelation of volatility, simulation B v. WIG data since 1991 
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Figure 3: The empirical (pluses) and simulational (crosses) autocorrelation func- 
tions of absolute log-returns vs. time. Autocorrelation functions obtained from 
simulations are consistent with the corresponding ones for market data |10j for 
both indicies. For comparison, simulations with no agents' esteem are also shown. 
The power-law relaxations of all autocorrelation functions are well seen over several 
decades. Detailed descriptions of the plots are in their titles and legends. 
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Log-returns histogram, simulation A v. WIG data since 1991 
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Figure 4: Histograms of daily log-returns both for market data [10] and obtained 
from simulations in the semi-logarithmic plots. The obtained results show the 
non-Gaussian properties of the price process, which is satisfactorily described by 
our model, in particular, for the stock market of a large capitalisation. Detailed 
descriptions of the plots are in their titles and legends. 
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